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particular, the pion form factor closely matches available experimental data. 
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1. Introduction 



Holographic models based on the gauge/gravity duality (descending from the AdS/CFT 
correspondence, originally introduced by Maldacena, see e.g., [Tp] ) have been studied and 
improved significantly over the last decade. Though still unable to provide a gravity dual for 
QCD proper, some of the signature features of QCD have been incorporated into semirealis- 
tic models. An important improvement was made within the model by Sakai and Sugimoto 
[3M] , who implemented a geometrical mechanism for chiral symmetry breaking into Witten's 
model [5] . The construction is based on the following strategy. The starting point is the near 
horizon limit of a geometry generated by a large number of color D-branes. To include 
fundamental quarks into the dual gauge theory one has to add a stack of Nf fiavor D-branes 
to this background. Strings stretching between the color and fiavor branes naturally trans- 
form in the fundamental representation U{Nc) x U{Nf). In case Nf -C Nc, the backreaction 
of the fiavor branes on the background can be safely neglected, i.e., one works in the probe 
brane approximation [6]. This corresponds to the so-called quenched approximation of (large 
Nc massless) lattice QCD, in which the backreaction of the fiavors on the colors is neglected 
completely while the color dynamics and its effect on the fiavors are accounted for. 
In order to have a geometric mechanism for chiral symmetry breaking one needs to add 
another stack of Nf anti-D-branes associated in the dual gauge theory to anti-quarks in the 
fundamental representation of U{Nf). This has to be accomplished in such a way that in the 
UV the two stacks are separated from each other, resulting in the dual gauge theory enjoying 
the full UilNf) X Uii{Nf) fiavor symmetry. In the IR however, if the two stacks of branes 
and anti-branes merge smoothly into a single stack, only a diagonal subgroup U£,{Nf) of the 
full fiavor group survives and chiral symmetry is broken. 

In the Sakai-Sugimoto model, this was beautifully realized by introducing Nf D8- and D8- 
branes into the near extremal geometry generated by a stack of Nc Z)4-branes in type IIA 
superstring theory. In order to obtain a 3 + 1 dimensional gauge theory in the IR, the X4- 
direction is compactified on a circle. The two-dimensional submanifold of the background in 
the (x4, r)-directions turns out to be cigar-shaped, leading to a U-shaped profile of the tra- 
jectories of the DS-DS-hiajies in these directions, thus constituting a geometric realization 
of chiral symmetry breaking as described above. 

The Sakai-Sugimoto model has been applied to a large number of different aspects in hadronic 
physics. For a review of mesons in this and other holographic models, see |^ and references 
therein. Form factors of vector and axial-vector mesons and pions in the Sakai-Sugimoto 
and other holographic models have been obtained in, e.g., PI^1U|11] . 

Although the Sakai-Sugimoto model was the first example of a stable geometric model of 
chiral symmetry breaking with completely broken supersymmetry, it suffers from shortcom- 
ings inherited from Witten's original model. Most notably, the model is only well-defined 
as a 3 -|- 1-dimensional field theory in the IR. As one goes to higher energies, KK modes 
become important at masses of the order of the inverse radius of the S^. At even higher 
energies, the dilaton (or string coupling) starts to blow up and one has to S-dualize. The UV 
completion of the theory is in fact the six-dimensional A^^^i (2,0) superconformal theory 
[T2IT3]. To overcome some of these problems, Kuperstein and Sonnenschein |T3] managed 
to incorporate stacks of fiavor D7- and D7-branes into a background with constant dilaton. 
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namely the Klebanov-Witten model [15], which is based on the singular conifold geometry. 
This embedding was engineered such that it breaks the remaining supersymmetry of the 
background completely and the two branches join smoothly in the IR|j 
In this paper we will focus on studying phenomenological aspects of the KS model [11], such 
as the spectrum of vector and axial-vector mesons, their mutual couplings and couplings to 
the pion as well as form factors of pions and (axial-) vector mesons. Our analysis indicates 
that vector meson dominance is realized in the KS model which means that the electromag- 
netic interaction of hadrons is mediated exclusively by vector mesons. The mass spectrum 
of (axial-) vector mesons was first obtained in [H], and our findings agree with the mass 
spectrum presented in that article. Many of the results presented below are quite similar to 
the results in the Sakai-Sugimoto D4 — D8 brane model [3P|8] . This is due to the similar 
shape of the Z)7-brane embedding. For the pion form factor, our results in fact reproduce 
the experimental data slightly better than the Sakai-Sugimoto model, as we will demon- 
strate below. Besides obvious numerical differences, one crucial difference is that, in the KS 
model, instead of a Kaluza-Klein mass scale Mkk, there is a mass scale M^, depending on 
the parameter tq, the minimal distance of the D7-branes to the singularity at r = 0, which 
also parametrizes a one-parameter family of solutions to the classical equations of motion of 
the D7-brane embedding. 
The paper is organized as follows: 

In section [2], we give a brief review of the Kuperstein-Sonnenschein model, followed by a 
discussion of vector and axial- vector mesons in section [3l Section H] contains a discussion of 
the pions and vector meson dominance. In section [5l we present our results for the various 
form factors. We formulate our conclusions in section [6l 



2. Review of the Kuperstein-Sonnenschein model 
2.1. D3-brane background 

The KS model is based on the DS-brane background with a conical singularity in type 
IIB superstring theory first studied by Klebanov and Witten [15]. This is an extension 
of the original AdS/CFT correspondence by studying strings on AdS^ x X^, where no 
longer is the 5-sphere but a five-dimensional Sasaki-Einstein space, namely the coset space 
T^'i = {SU{2) X SU{2))/U{1) with Einstein metric 

1 / ' \' 1 ' 

\ i=l / t=l 

Before taking the near-horizon limit, the six dimensional internal space transverse to the 
Z?3-branes is a Calabi-Yau manifold Yg with a singularity, called a conifold, which is a cone 
over the base manifold T^'^ (topologically, T^'^ is 5*^ x S'^\ Namely, 

c^S(io) = -^Vf^udx^dx" + -^ds\y (2.2) 

lot of work has been done studying supersymmetric embcddings into the Klebanov-Witten and 
Klebanov-Strassler models (e.g., }16|17|18] ). However, it turns out that these embeddings only allow a 
single branch of branes and thus they are not useful as models of chiral symmetry beaking. 
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where 

o?4) = (dr^ + r^rfsli.i) , (2.3) 
and the AdS^ radius is = ^ngsNJ'^. 

The field theory dual resulting from a stack of A''^ il'3-branes placed at the conifold singularity 
is a A/" = 1 superconformal field theory with gauge group SU{Nc) x SU{Nc) and global 
symmetry {SU{2) x SU{2)) x f/(l). It contains four chiral superfields, Ai,i = 1, 2, a doublet 
of the first SU{2) factor transforming as {Nc,Nc) and Bj,j = 1,2, a doublet of the second 
SU{2) factor transforming as {Nc,Nc). The theory has a marginal superpotential W = 
e'^e^'^AiBjAiBm, and the R-charge of all chiral superfields is |. 



2.2. D7-brane profiles 



One of the achievements of [H] was to study a DT-brane configuration spanning the space- 
time coordinates x^, the radial coordinate r and the parametrized by the one- forms /j. 
The one-forms /j, i = 1,2,3 are defined as (for details, cf. 



(2.4) 





where the SU{2) left-invariant Maurer-Cartan forms w'^ji 
following condition: 



. 3 

-y 

2 ^ 



1,2,3 are determined by the 



(2.5) 



The explicit form of the one- forms w'^ and the definition of the matrix X can be found, e.g., 
in [36] • It should be noted that the structure involving the 5*0(3) matrices determines the 
fibration of the over the S*^. This leaves two transversal coordinates, namely those on the 
S'^, 9 and (f). We will assume that 9 and do not depend on the coordinates (note that 
our ansatz preserves one of the SU{2) factors of the global conifold symmetry). This can 
be justified by expanding the action around the solution and verifying that contributions 
from non-trivial modes only appear at second order in the fiuctuations. Therefore the 
classical solution will only depend on the radial coordinate, 9 = 9{r), (p = (f){r). Plugging 
this ansatz into the DBI action for D7-branes yields the following Lagrangian (for details, 
see the discussion in [S]): 



r 



— + sm2 9 
or I \ or 



(2.6) 



We can use the SU{2) invariance of the Lagrangian to set 6* = 7r/2, i.e., restrict the motion 
to the equator of the S"^. Then the solution of the field equations is given by 



cos I 



rox4 
r 



(2.7) 
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The trajectory 0(r) has two branches, namely G [— 7r/2, 0] and G [0, 7r/2]. In the extremal 
case ro = 0, there are two r-independent solutions, (j)± = ±^7r. For tq ^ 0, the solution 
(j){r) starts at 0(r = tq) = and approaches the asymptotic values 0± as r — )■ oo. 
In conclusion, [13] found a one-parameter family of solutions parametrized by ro sharing the 
same boundary values (f)±. For tq > the conformal symmetry of the underlying Klebanov- 
Witten background is broken, since there is no AdS^ factor in the induced metric. It is 
worth noting, however, that in the limit r — )■ oo, one recovers the AdS^ metric and conformal 
invariance is restored. As a consequence of this asymptotic behavior, we expect that the KS 
embedding reproduces QCD results better than the Sakai-Sugimoto model in the UV regime. 
Introducing the coordinates y = r"^ cos (^'^'Pj z = r'^ sin ^^^d noting that y = Tq is 

constant along the classical configuration solving (12. 7p . the induced metric on the D7-branes 
can be written as 



^^(8) = ^,%.dx^dx^ + \^-^^dz- - y^dzh + -Ul + /I) + Y j ' (2-8) 
where = + Tq determines the relation r = r{z). 
3. Vector and axial- vector mesons 

The strategy for calculating vector meson spectra utilized here will closely follow |8]. Ac- 
cording to the gauge/gravity correpondence, mesons of the dual gauge theory arise as modes 
of open strings stretching between the D7- and Z)7-branes. In particular, vector and axial- 
vector mesons correspond to fluctuations of the U{Nf) gauge fields living on the Z)7-branes. 
To first order, these fluctuations are captured by the action 

5-^7 = -(TraOVr j d'^xdzd^Qe-" y^^M^MNF^'"") + Scs, (3.1) 

where the first term is the Maxwell term of the non-abelian Dirac-Born-Infeld actioij^ of the 
D7 branes and the second term is the Chern-Simons term, which is irrelevant for the topic at 
hand. The non-abelian field strength Fmn = QmAn — OnAm + [Am, An] contains the gauge 
field components Am = {A^, A^, Aa), /i = 0, 1, 2, 3, a = 5, 6, 7 and g(^s) is the determinant of 
the induced metric on the -D7-branes, given by (12. 8p . In what follows, we will set Aa = 
and assume that the other gauge field components are independent of the coordinates. 

Five-dimensional effective action. Integrating the D7-brane action (13.11) over the 
coordinates leads to a five-dimensional effective action for the fiuctuations of the gauge 
fields, 



S5d,eS = -K I d X I dz tT 



(3.2) 



^We neglect higher order contributions to the DBI action which are suppressed by higher powers of a' . 
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where n = 2(7ra')^/i7(7^ ^^VolS^ = Here we have defined the dimensionless variable 



z = ^ and an ro-dependent mass scale M^(ro) = The functions 



'0 



C{z) = {z^ + 1)-^/' and D{1) = {l^ + 1)^/^ (3.3) 



are rescaled, dimensionless variants of the functions Ciz) and Di^z) introduced in [T] 
eq. (4.14). It turns out to be convenient to work in the = gauge in which can 
be expanded in the following way: 

oo oo 

A^(x,5) = V^,{x)^A^{x)Mz) + ^t/")(x)V^2n-i(5) + 5^a(r)(x)V^2n(5), (3.4) 

n=l 71=1 

wher^ ^0(5) = ^r'li/l) ^ 2-^1 (|j f, f, —5^), and the fields v]f^ and a^""* represent vector and 
axial-vector mesons, respectively. The pion field n(x) appears in the expansion of 

v.(^) = ^(4;(^) + 4-.(^)) (3-5) 

= \u {^lM + Q,) it + \i- (^h,m(x) + a^) iz\ 

M^) = ^ (4;(^) - 4.(^)) (3.6) 

= (^L,(X) + 9,) - + 9,) c\ 



where Al^{x) and y4/j^(a:) are external gauge fields and we have defined = ^_ := e* 
The normalization conditions and equations of motion for the wave functions ipn{z) can be 
written as: 



K J dzC{z)tjjm{z)tlJn{z) = 5m,n, (3.7) 

-{C{z))-^d,{D{z)d,M^)) = XnM^), (3.8) 
with A„ := M^-g, = M^/Ml 



^The normalization constant is obtained by demanding 



Jim A^(x^i) = 4;(a.^), Jim A^{x^^,~z) = A%- {xn- 
In the Az ~ ^ gauge, we can write 

00 00 

A,{X^\~Z) = 4;(x^)^+(f) + A%-^{xni^.{~z) + ^ <)(x)^2„-l(5) + ai"\x)^/j2n{z), 

n—1 n—1 

with V'±(^) = 5(1 i ^/;o(5))- The normalization constant now follows from imposing 



lim 4'o{z) = ±1- 

i— >±oo 



Note that dzip± = zt^dztpo{z) oc D{z)' 
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Four-dimensional effective action for the mesons. Substituting the gauge field (13. 4p 
into the action S^d,cs^ eq. (13. 2p . we find the four-dimensional effective Lagrangian for the 
vector and axial-vector mesons (ignoring divergent terms from non-renormalizable contribu- 
tions), CidfiS = ^Ad,div + ^^=2^4d,cS- ^hc kiuctic part of the Lagrangian reads 

42efr = ^tr(9,e^-94"))2 + ^tr(9,4")-94"))2 + tr(z9,n + A^,f 

+M.^.tr (") - ^V,) ' + Mf„tr - ^ a) ' , (3-9) 

where we have used the following redefinitions of meson fields in order to diagonalize the 
kinetic terms: 

Moreover, we define 

= X2n^lMl Ml„=\2nMl (3.11) 

V, = ]^{Al, + Ar,), A, = \{Al,-Ar,). (3.12) 

The constants 

/ + 00 
dzC{z)ll)2n-l{z), 
-oo 

= -2kM^ {D{z)d,^lj2n-i{z)) , (3.13) 

r+oo 

ga« = kM^„ / dzC{z)i;2n{z)M^), (3-14) 



are the couplings between a massive vector meson wf,"^ (axial- vector meson ai"^) and an 
external U{1) field representing a photon (an external axial U{1) field A^). 
The interaction part of the effective Lagrangian governing the interactions among vector and 
axial- vector mesons is given by 

+ g,.„{d'^d^'^> - d''d^"^n H'\ 4"^^] - S^r^]) } . (3.15) 
The 3-meson coupling constants are given by 

g^e^m^n = dzC{z)lp2£-l{z)'ip2m-l{z)lp2n-l{z), (3.16) 

g^la^a" = 1^ dzC{z)4j2e-l{z)tp2m{z)^p2n{z). (3.17) 



The next step is to numerically calculate the wave functions ipn{z) in order to be able to 
obtain the couplings gyn^ gytym^n, g^ia^a" ^^id the masses and M2„. 
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3.1. Numerical study of the wave functions 

To study and numerically solve the field equations for the vector and axial-vector me- 
son modes, we will employ a shooting method (cf. [31I8] ). In order to find the correct 
asymptotic behavior of the wave functions ipn{z), we expand them into a Frobenius series 
ipniS) = z~'^Ylim,^n,m.z~'^ for large z. Substituting this expansion into (I3.8p . and imposing 
the normalization condition (13.71) . we find that asymptotically il)n{z — )■ -|-oo) ~ 
Therefore, it will be convenient to define ipn{z) = z^^'^^niz). Introducing the new variable 
5 = e'', the equation of motion for ipn can be cast into the following form: 



where 



d^^^Pn + Ad^tPn + Bi^n = 0, (3.18) 



A{~z = e^) = -^^^^ = E^^^"'^'' (3-19) 

q ~-2 ~-l/2 °° 

The first few coefficients are 

Ao = -^,Ai = A2 = A3 = 0, ^4 = As = = ^7 = 0, As = I . . . (3.21) 

3 5 3 

Bo = 0,Bi = A„, i?2 = -63 = 0, B4 = -, i?5 = — ^A„, Bq = B7 = 0, Bg = — -, (3.22) 

Finally, by expanding 

00 

^„(5) = ^a^5-™/2, (3.23) 

m=0 

with ao = 1; and substituting this expansion into f l3.18p . we find the following recursion 
relation for am'- 

ai = -2Bi, (3.24) 

-1 / T m—l m—l 



^ + ^ j I 2 X/ kAm-kOCk - 2_, Bm-kOCk 1 • (3.25) 

^ \ k=l k=0 ' 



These relations can now be taken as input data to solve the field equations (13. 8p numerically 
via a shooting method. 



3.2. Results 

Parity. Let us briefly discuss the behavior of the wave functions under the parity transfor- 
mation (t, X, z) —J- (t, —X, —z). A parity transformation exchanges the chiral components of 
the external gauge fields, Ai f-)- Ar. The vector and axial-vector meson modes exhibit the 
following behavior under parity: 

t;(")(t, -x) = +v'-;\t, -x), a^;:\t, -x) = -a^;\t, -x). (3.26) 
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To ensure gauge invariance of the five-dimenional gauge fields, we impose the conditions 

'4'2n{-z) = +t^2n(5), 1p2n-l{-z) = -^2n-l(5)- (3.27) 

If we demand regularity of the wave functions ipn{z) at the origin z = 0, from the above 
conditions we conclude that 

9i^2„(0) = 0, V^2„-i(0) = 0. (3.28) 

Using the large z behavior, eq. (13.231) . and the conditions at 2 = 0, we are now able to 
numerically solve the equations of motion (13. 8p . We found solutions ipn{z) and the corre- 
sponding eigenvalues A„ for n = 1, . . . , 60. The results for a selection of wave functions are 
shown in figure [TJ For large z one finds the expected behavior. In the small z limit. 
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Figure 1: Wave functions y/Kipn{z) for n = 


= 1,2,3,4,30,59. 



z 



the equations of motion reduce to 

<(5) + A„7A„(5) =0, (3.29) 
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which has sinusoidal solutions 



'tp2n ~ sin(v^A^2;), ilJ2n-i ~ cos(v/A2n-i5) (3.30) 



showing oscillatory behavior. This is compatible with the observed behavior of ipn{^) in 
the small z regime. The corresponding mass eigenvalues are summarized in table ^ The 
numerical values found here are consistent with the values quoted in [M], eq. (4.24) jj 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 




0.258 


1.384 


3.432 


6.393 


10.266 


15.049 


20.745 


27.348 


34.866 


A2n - A/2 


0.691 


2.292 


4.798 


8.215 


12.408 


17.783 


23.932 


30.995 


38.920 



Table 1: Some numerical values for the dimensionless vector and axial- vector masses. 

The lightest meson, v^^\ is a vector meson that could, for the sake of comparison with QCD, 
be identified with the p-meson p(770). The second lightest meson, a^^\ is an axial- vector 
meson identifiable with ai(1260), the third lightest, vector meson identifiable with 

p(1450) and so on. In table [21 we compare the mass ratios of the first few mesons in the KS 
model to the Sakai-Sugimoto (SS) model and experiments (QCD). It is important to keep 
in mind, however, that the predictions of any gauge-gravity model are stricly speaking only 
valid in the large Nc, large t'Hooft coupling limit and that we are not studying a QCD-dual 
here, but rather a dual to the quiver gauge theory described above. 



k 


( ^k + l \ 

V ^1 Jks 


V ^wss 


\ ^1 


) QCD 


1 


2.68 


2.4 




2.51 


2 


5.36 


4.3 




3.56 



Table 2: Ratios M|_^^/M^ = A^+i/Ai for = 1,2. M\ corresponds to the mass squared of 
the p-meson, which can be matched in our model by setting Ml = (1527MeV)^. This can 
be done by adjusting appropriately. 

Regge trajectories. Next, we examine the dependence of the eigenvalues A„ on the radial 
number n, which we will refer to as a Regge trajectory. These eigenvalues determine the 
meson masses according to = A„M^. The Regge trajectory for both the vector and 
axial- vector mesons is shown in figure [21 We also display a logarithmic plot with linear best 
fits for both small n = 1, . . . , 5 and large n = 6, . . . , 60. In the small n region, the best 
fit is —1.41 + 1.61 logn, corresponding to A„ ~ ra^'®^, while in the large n region, we find 
—1.97+1.95 logra, corresponding to A„ ~ n^-^^. This result is consistent with previous results 
presented in the literature, where the Regge trajectory for vector and axial- vector mesons 

"'in order to compare the numerical values, it is important to observe that the A„ of [14| in reality are 
differently normalized meson masses For direct comparison between the two papers, one has to use the 

relation A„,horc = -yf^- 
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Figure 2: Regge trajectory for vector and axial vector mesons. In the left graph, A„ is 
plotted versus n, for both vector mesons (+) and axial vector mesons (-). The right graph is 
a logarithmic plot of both vector and axial- vector mesons. The straight lines are two linear 
best fits for small n < 5 and n = 6, . . . 60. 



was found to exhibit asymptotic quadratic behavior A.„ ~ in the DA — DS model [Sj , the 
hard- wall model [T^2U|21|22|23] and the D3 — D7 model Thus, this seems to be a 

general feature of holographic models of mesons. 

Decay and coupling constants. We are now ready to calculate the decay and coupling 
constants. Since we will be interested in studying (elastic and non-elastic) form factors, i.e. 
interactions of particles with external photons, which naturally couple to vector mesons, 
we will focus here on the vector meson decay constant g^n and the coupling constants be- 
tween the intermediate vector meson and the external vector mesons g^n^i^m,, and between 
the intermediate vector meson and the external axial- vector mesons g^n^i^m, respectively. 
The calculation was carried out using eqs. f l3.13p . fl3.16p and fl3.17l) . as well as the numer- 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 




0.829 


2.965 


5.864 


9.347 


13.331 


17.759 


22.591 


27.793 


33.345 




0.383 


-0.110 


0.00176 


-0.000210 


0.0000631 


-0.0000194 


0.0000265 


0.0000311 


-0.0000213 




-0.110 


0.258 


-0.119 


0.00254 


-0.000349 


0.0000621 


-0.0000324 


-0.0000320 


0.0000225 




0.258 


0.0579 


0.136 


-0.131 


0.00382 


-0.000526 


0.000108 


0.0000210 


-0.0000217 




1.783 


4.337 


7.539 


11.281 


15.294 


20.127 


25.147 


30.529 


36.243 




0.270 


0.156 


-0.122 


0.00351 


-0.0000443 


0.0000673 


-0.0000103 


-0.00000183 


0.00000236 


y/^9v"a^a'2 


-0.115 


0.147 


0.138 


-0.128 


0.00384 


-0.000507 


0.0000805 


-0.00000579 


-0.00000342 


y/^9v"a'^a'^ 


0.249 


0.0540 


0.0484 


0.117 


-0.136 


0.00422 


-0.000594 


0.0000829 


0.00000500 



Table 3: Dimensionless decay and coupling constants for vector and axial-vector mesons. 

ical results for the normalized wave functions found above, for n = 1, . . . ,60. The results 
of the computations are shown in table [3] for both vector and axial-vector mesons. For 
completeness, we also include the coupling constants ga"- 
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4. Pions 



Skyrme model. Let us briefly discuss the effective Lagrangian for tlie pions coming from 
tlie DBI part of tlie D7-brane action. To tliis end, it is convenient to omit tlie meson fields 

vjr\ alj^\ and to work in a gauge where = 0, ^.(x'^) = 1 and ^.^^{x'^) = U{x^) := . 
Then we are left with the mode expansion 

~z) = U-\x'')d,U{xn4>+{z), (4.1) 

which leads to 

F,, = [U~'d^U,U~'dM]i^+{^+-l), (4.2) 

F,, = U-^d,Ud,^+. (4.3) 

The resulting action for the pions can then be written as 

Sfd^K = ^ j d'x j dzii{^-C{~z),ljl{^+-lf[U-'d,U,U~'d,uf + MlD{z){^^^ 

(4.4) 

This is precisely the action of the Skyrme model 



^Skyrme = j d^X (^^tl {U-'d^Uf + ^tr [U-'d,U, U-^d,U\ 



(4.5) 



with the following identifications: The pion decay constant /tt (~ 93 MeV experimentally in 
QCD) is given by 

fl = ^^^Ml j dzD{~z){d,^^f = -^^^^^M^ « 0.76/.M,^ (4.6) 

and the dimensionless parameter is calculated as 

ef = 32k j dz]^C{z)^l{ip+ - if = k j diC{z){l - i^lf ^ 3.93k. (4.7) 

Here we used the definition of ip+{z) = |(1 + ipo{^))- In order to match the experimental 
value for f^,, the dimensionless constant k has to be adjusted accordingly: 

K ^ 4.87 X 10"^ (4.8) 

In large Nc QCD, one expects ~ 0{Nc) and e^^ ~ 0{Nc), which is consistent with the 
dependence on k (cf. the expression below eq. (13.21) ). 

Interaction terms involving pions. So far, we have only studied the part of the interac- 
tion Lagrangian C\J^q involving mesons. For the computation of the pion form factor below, 
we also need to introduce the part involving pions and mesons simultaneously. To this end, 
we rewrite the kinetic term for the vector mesons using the redefinition 

t;J^)=4-) + ^[n,9,n]. (4.9) 
12 



Now, written in a form closely resembling the effective Lagrangian given in the literature 
(cf. [IpG] ). the relevant terms are 



+M^.tr(^;)-^V.)%M:.tr(.H-^^.)^ 

-2^....tr(t;^)[n,a'^n]) + (2^^:|f^ - 2^ tr(v^[n,9''n]) 



+^tr[n, d,ur - 2^tr[a,n, d.Uf 



Here we have introduced the coupling constants 



71^2 r+Qo 
K-^ / rf5C(z)^2„-l(5)(l-^o') 



d~z^o{~z)d,^2n-i{~z). (4.10) 



Vector meson dominance. Using the equation of motion fl3.8p . together with the com- 
pleteness relation 

oo 

nY,C{i')Mz)M^')=^{~^-~z')^ (4.11) 

n=l 

one can prove the important sum ru 1^ 

oo 

E^^-^^ = l- (4-12) 

n=l 

This shows that the Vvftt coupling in the above Lagrangian vanishes, leading to vector 
meson dominance in the pion form factor. This result descends from a more general class of 
geometries and brane embeddings leading to a five dimensional effective action as in (13. 2p . 
General C{z) and D{z) terms need to obey a few conditions in order to reproduce vector 
meson dominance and other properties of pion and meson phenomenology. A finite pion 
decay constant requires that (cf. f l4.6p ) 



/. ~ j dl— (4,13) 

be finite. Renormalizable five dimensional gauge fields require vanishing wave functions at 
z — )■ ±oo whose behavior depends on C{z) and D{z) in (13. 7p . Parity requires that both C{z) 
and D{z) be even functions of z. Noticing that the equations of motion (13. 7p are of Sturm- 
Liouville type, the above requirements set the conditions for an orthonormal set of wave 



^This sum rule turns out to be universal in holographic models. It has be shown to hold in the soft-wall 
and hard-wall models [T^, USZ-DZ-models [5] and the Sakai-Sugimoto Z34/Z)8-model [3]. 
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functions ipn{z)- Vector meson dominance hence follows from applying the wave functions 
completeness relation (14. lip and equation of motion (13. 8p to the sum rule expressions for the 
pions ( I4.12P and mesons (IS.lOp . In [9], the authors claim that vector meson dominance is a 
feature of asymptotic AdS spaces, but it should be noted that it is in fact more general and 
also applies to other geometries, e.g., the Sakai-Sugimoto model. 
Similar to (14.121) . one can prove two more sum rules: 



oo 2 

9v"7 



E 



n=l 



1 



1 1 



oo 2 

^ "4e|/f 



(4.14) 



In conclusion, the kinetic and interaction Lagrangian involving both mesons and pions can 
now be written in its final form: 



(2) 



4(i,efr 



C 



(3),meson+n 
4d,eff 



itr(9,e^ - d.v(;yf + ^tr(9,aW - d^d'^^^y 
+tr(z9,n + UA,f + Mlii (^trf _ ^v,) 

+M2„tr - ^ a) - 2^?.n..tr(i;^) [H, d^U]) 



3/; 



,tr[n,a^n]2- 



tr[a^n, d,Ii] 



(4.15) 



The wave functions calculated in section [3?T] fall off rather slowly at large z. This results in 
fairly slow convergence of the relevant integrals involved in studying the pion form factor. 
For the results presented below, we have integrated numerically up to 5max = 10^^ using 
Mathematica. With this precision, we can roughly trust the first 15 coupling constants 
dv'^nn, n = 1, . . . , 15 introduced in eq. (I4.10p . In table HJ we list some numerical results for 
these coupling constants. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 




0.335 


-0.0425 


0.00914 


-0.00216 


0.000545 


-0.000141 


0.0000535 


0.0000155 


-0.0000140 



Table 4: Coupling constants for a vector meson t>" coupling to pions. 



5. Form factors 
5.1. Pion form factor 

The electromagnetic form factor represents the interaction of a particle with an external 
photon. In this section, we will study the pion form factor F.,^{q'^) (cf. figure 3), which is 
defined by the following matrix element of the electromagnetic current, 

{7r'^{p)\J^m\Ap')) = r'\v + VrmP-v'?). (5.1) 
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Figure 3: Feynman diagram for the pion form factor, show- 
ing vector meson dominance. 



where is the conserved vector current coupled to and /"'"^ are the structure constants 
of U{Nf). By combining the relevant vertices and propagators (cf. eq. fl4.15p ). according to 
figure 3, we arrive at the following expression for the pion form factor, 

oo 

^-tf) = E|Ti^- (5-2) 

Comparing eq. (14.121) to eq. (15. 2p . one finds that Ft^{{]) = 1. We can use this as a test of the 
quality of our numerics. Summing up the first 15 terms in (15. 2p . we find 
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- j^9v-nn ~ 0.999902. (5.3) 



n=l 



The numerical results for the pion form factor are presented in figure HI 

To study the large behavior of the form factor, we perform an expansion in powers of q~^: 

^ n=l V 9' / 

We can calculate the coefficients of the leading term, using our numerical results: 

15 

^ g.ng.r.^^ ^ 0.443(GeV)2, (5.5) 

n=l 

which is (to a high accuracy) the asymptotic value of x F^^lq^). The large behavior 
of F^^lq"^) ~ is a reflection of the fact that pions (as mesons) are made up of two 
constituents (see the discussion in [11], appendix D, specifically eq. (Dll)). Note that the 
large g^ behavior corresponds to — >■ 0, i.e., taking the conformal limit ro — >■ in the KS 
model. It is also interesting to observe that for large g^, the slope of g^ x -F7r(g^) in figure 
m depends on ^n9v"9v"mTM^'^, which converges to zero quickly in the KS case, signaling 
conformal behavior, while in the Sakai-Sugimoto model it turns out to be nonzero and 
negative (~ — 2.5(GeV)^). Its slope will still eventually converge to zero for very large g^ 
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Figure 4: The plot on the left shows the pion form factor for the KS model (solid blue 
line). On the right, we present a plot of x F^.i^q'^), which exhibits the expected 
behavior at large g^. The KS model (solid blue line) agrees very well with the data obtained 
from QCD experiments, taken from [271128] . The error bars represent all experimental and 
analysis errors, but not "model uncertainty", as explained in the aforementioned papers. 
For comparison, we also included the results for the Sakai-Sugimoto model (solid red lines), 
compiled using data from [HIS] . Note that the results for the Sakai-Sugimoto model should 
only be trusted up to the KK scale M|^^ ~ (0.946GeV)^. The dashed red lines represent 
extrapolations based on our numerical calculations in the Sakai-Sugimoto model. 



since it comes with a factor of q ^. 

We can also numerically check the validity of the sum rules fl4.14p . 

15 2 15 2 

.2 _ 0.333285, elf^ V ^ 0.249927. (5.6) 

n=l " n=l " 

Moreover, we are now in a position to estimate the charge radius of the pion. Expanding 
F^^^q^) for small q^ as 

oo 
n=l 

we compute 

n=l 

Summing up the first 15 terms, as before, yields 

{rl) ^ 24.707M-2 = (0.641fm)^ (5.9) 
which is close to the experimental value of (r^)exp ~ (0.672fm)^. 



5.2. Vector and axial- vector meson form factors 

Here we will calculate the form factors for vector and axial-vector mesons, again following 
the strategy outlined in [8] for the D4-D8 model. Using a sum rule similar to eq. fl4.12p in 
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the case of pions above, 

oo 

^ ] 'T72~9v'^v"^v' = ^m£, (5.10) 
n=l 

one can show that all photon-meson-meson couplings are cancelled in our model. Because of 
this the photon interacts with a meson only through intermediate vector mesons, which is a 
realization of vector meson dominance (VMD) in electromagnetic scattering in this model. 
Numerically, the convergence of the relevant integrals is even slower than in the pion case. 
However, we find a vector meson dominance pattern that is very similar to the one described 
in [To] in the AdS/QCD hard wall model, where the dominant contributions to F^i{q^ = 0) 
were shown to come from the first two bound states (to 10"'^ accuracy). Here we find, e.g, 

6 

F,^{q^ = 0) ^ T^9v^v^^^ = 0.999863, (5.11) 

n=l 
6 

F„2(g2 = o) ^ ^ = 0.999840. (5.12) 

n=l 

In order to compute the various form factors we first obtain a general expression valid 
for the elastic and non-elastic cases. Then, for the elastic case, we calculate the electric, 
magnetic and quadrupole form factors and briefly touch upon the large behavior of the 
longitudinal, transverse and longitudinal-transverse form factors. Also, for the first vector 
and axial- vector excitations, p(770) and ai(1260), we compute the electric radius and the 
magnetic and quadrupole moments. We conclude this section with a brief discussion of the 
so-called transition or non-elastic form factors. 




Figure 5: Feynman diagram for vector meson form factor. 
A similar diagram holds for the axial- vector meson replac- 
ing the external lines by a^"*-*, a^^-*. This graph 
also exhibits vector meson dominance. 

As before in the pion case, the form factors are calculated from the matrix elements of 
the electromagnetic current. The interaction of a vector meson with an off-shell photon is 
given by the matrix element 

{v^^'^\p),e\JAq)\v^'^\p')A) = (27r)V(p'-p-g)(i;("^)'^(p),e|J^<^(0)|^;W^(p'),e'), (5.13) 
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where v^'^^ and v^^^ are the initial and final vector meson states with momenta p and p' = p+q 
and polarizations e and e'. The operator J'^ is the Fourier transform of the electromagnetic 
current J'^{x). This matrix element can be calculated from the corresponding Feynman 
diagram shown in figure 5. From the effective Lagrangian f l3.9p . together with the interaction 
terms fl3.15p . we find 

e| J^^(0)|^;W^(p'), e') = e'^e'T'^ [vUq - p)p + Vup{2p + q)^ - Vpaip + 2g),] 



oo 



g2 + M2„ 



(5.14) 



n=l 

where M„n is the mass of the vector meson w*^"). Defining the generalized vector meson form 
factor as 

oo 

Fw(g) = )i^^^^^, (5.15) 

n=l " 

and using the sum rule eq. (15.101) . we find 

(^('")'^(p), 6| J^^(0)|t;W''(p'), e') = e^d'r"' ^.(g - p), + r/,,(2p + g). - r/,.(p + 2q\\ 

X { (^'^'^ - ^) ^.-.^(?') + '^-^^} • (5-16) 

Taking into account the transversality of the vector meson polarizations, e ■ p = = e' ■ p', 
and noting that the term involving the factor bmi does not contribute since, in the elastic 
case m = £, we have 2p ■ g + = 0, we obtain 

(?;("*)'^(p),e|J^^(0)|i;(^)'(p'),e') 

= [^.p(2p + g). + 2(r/.,g, - r^^^g,)] (^r/^'^ - F,^..(g') . (5.17) 

Note also that this matrix element satisfies the tranversality condition g^(i;'^™'^|J^(0)|w''^^) = 0. 
For axial-vector mesons a similar calculation can be done to find the form factors from the 
matrix element (a("*)"(p), e| J^"(0)|a(^)^(p'), e')- In this case we should replace the external 
vector meson lines by the axial- vector mesons a*^"*^ and a*^^-* in the Feynman diagram of figure 
5. The internal vector meson line f remains unchanged due to vector meson dominance. 
This is ensured by the sum rule (compare to eq. fl5.10p ): 



oo 



-TlPrQv-^a^at — ^ml: (5.18) 

n=l 

Again, we can test this sum rule numerically, adding up the first 6 terms of the corresponding 
sums (for a*-^^ we sum up to 7 in order to improve numerical accuracy), 

6 

i^ai(g' = 0) ^ 5^ ^Wai = 0.999872, (5.19) 

n=l 
7 

F,2(g2 = 0) ^ ^ = 0.999804. (5.20) 

n=l 
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Thus, the generahzed axial-vector meson form factor is 



(5.21) 



Elastic case. In order to obtain the elastic form factor for vector mesons, we consider the 
previous calculation imposing that the initial and final vector meson states v^"^^ coincide. 
Then, from eq. (15.170 . we find 



= r'^iie ■ e'){2p + qr + 2 [e^{e' ■ q) - e'^e ■ q)] ] F^^{q') 



where F^m^q^) is the elastic form factor: 



F^m {q^ 



n=l 



(5.22) 



(5.23) 



Similar relations hold for the axial- vector mesons, merely replacing v^"^^ by a^'^\ We calcu- 
lated numerically these sums from n = 1 to n = 6 using the results obtained for the masses 
and couplings. The mass scale is = 2.332 (GeV)^, as explained above in section |3l We 
plot in figure E] the elastic form factors for the vector meson p(770) (v^^'^) and axial-vector 
meson ai(1260) (a^^^). Note that when — )■ 0, the vector and axial- vector form factors go 
to one, due to the sum rules (15.101) and (I5.18p . respectively. 




8 10 12 14 



q'(Gey-) 



Figure 6: Elastic form factors for for vector mesons (left) and axial-vector mesons (right) for 
the first two excited states: i = 1 (solid line), i = 2 (dashed line). 

It is interesting to investigate the large behaviour of these elastic form factors. Per- 
forming an expansion in powers of g~^: 



1 °° / 

i^^m(g^) = — ^ ^^ gv^gv^v^y-^ ( 1 



n=l 

oo 



F,™(g2 



n=l 



ML 



f + 0{q~''] 
+ 0{q~'] 



(5.24) 
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we see that the dominant terms would be of order q ^. We calculated the coefficients of 
these terms using our numerical results with n = 1, 6 (or 7, respectively). We found 

6 6 

Y,9v^9v-v^v^ ~ 0.000176(GeV)2, ^ ^ 0.000220(GeV)^ 

n=l n=l 
6 7 

^^^n^^„„2„2 ^ -0.00309(GeV)^ ^^„n^„„„2„2 ^ -0.00546(GeV)^ (5.25) 

n=l n=l 

These results indicate that the super convergence relations 

oo oo 

^ ^ gv^gyn-y-^y-^ = , ^ ^ 9 9 a"^ a"^ = 0, (5.26) 
n=l n=l 

hold in the KS model. Thus, from eqs. fl5.24p we expect that the form factors decrease 
approximately as for large g^. In figure [TJ we plot the elastic form factors for the ffist 
and second vector and axial- vector states multiplied by q^. 
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Figure 7: times the elastic form factors for vector mesons (left) and axial-vector mesons 
(right) for the ffist two excited states: i = 1 (solid line), i = 2 (dashed line). 

This corroborates the expected fall-off with g""^. The deviation from the g""^ dependence is 
related to the non- vanishing of the numerical sums fl5.25p . 

Electric, magnetic and quadrupole form factors. It is instructive to calculate the 
electric (-Fe), magnetic (-Fm) and quadrupole (Fq) form factors from the previous results. 
These form factors can be defined as (cf. |10j ) 



2 



Fe — Fi + 



2 



Fm = Fi + F2, Fo = -F2+ (l-^)f3, (5,27) 



where Fi , F2 and F^ are related to the matrix element of the electromagnetic current for a 
spin one particle, 

(p,e|J^(0)b',e') = (e.e')i2p + qrF^iq')+[e^ie'■q)-e''ie■q)] [F^iq') + F,{q')] 

+ \{q-e'){q-e){2p + qrFs{q'). (5.28) 
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Then, from eqs. ( I5.22p and fl5.28p we find that for a vector meson v 



(m) 



= i^r^ = , = , (5.29) 

where F^m is given by eq. f l5.23p . Therefore the electric, magnetic and quadrupole form 
factors predicted by the KS model for vector mesons are 

Fe ^ = (1 + ^2)-^^™ ' ~ , Fq "* = — F„m . (5.30) 

The same formal results hold for the axial- vector mesons a^"'\ In the following, we will 
estimate three important physical quantities from these form factors which are associated 
with the vector mesons: the electric radius, the magnetic and quadrupole moments. 
Electric radius. The electric radius for the vector and axial- vector mesons are given by 

{rU = -6^i^r^(^^)l.^=o , (rl^) = -6Ai.(f'")(,2)|^,^„. (5.31) 

Using our numerical results for the form factors for the lowest excited states p and ai, we 
find the electric radii (setting = and = M^^, respectively): 



00 



n=l P 

00 



^ ML Ml ^ ' 



Our results are slightly smaller than the results for the electric radii obtained in the DA-D8 
model [8] and the hard- wall model [10]. 

Magnetic and quadrupole moments. Owing to the fact that Fyn and Fa^ go to one 
as — 0, the magnetic and quadrupole moments in the KS model formally reproduce the 
results obtained in other holographic models, e.g., the Sakai-Sugimoto model [S] or the hard 
wall model [TO] : 

/i := FM(g^)|g2=0, ^/i„n/„n=2 (5.34) 

1 11 

D ■■= — ^FQ{q^)\g2=o, ^ Dy„. = = (5.35) 

Of course, the numerical values of these quantities in the KS model differ slightly from 
other holographic models, since the masses Mj,»i and M^n differ from model to model. 
Decomposition in terms of transverse and longitudinal polarizations. It is in- 
structive to study the form factors of vector mesons with specific polarizations. The def- 
inition of the transverse and longitudinal polarizations and the corresponding form fac- 
tors can be found in, e.g., [2^IIB] . The large behavior of these form factors in the KS 
model is determined from F^™, eq. f l5.23p . We obtain the following asymptotic behaviors: 
F^^P ~ g-4 , F^^'P ~ ^ p^vp _ ^ agreement with the DA-D8 model [8] and QCD 
calculations (see, e.g., [29]). 
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Non-elastic case. Finally, we consider the form factors for the case that the initial and final 
meson states are different. These are called generalized, transition or non-elastic form fac- 
tors. The interaction with the photon is determined by the eqs. (15.15^ . ( I5.2ip . with m ^ i. 
We calculate these form factors explicitly for the initial vector meson state v^^^ and final 
states w*-'"^ with m = 2, 3, 4 and the same for axial-vector mesons. 

The results are plotted in figure [H compared with the corresponding elastic form factors. 
Here, denotes the momentum transfer. 




Figure 8: Transition form factors for vector mesons F^\^i (left graphs) and axial-vector 
mesons -Faia* (right graphs), for i = \ (elastic case, solid line), i = 2 (dashed line), 2 = 3 
(dot-dashed line) and i = 4 (dotted line). 
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Figure 9: times the transition form factors for vector mesons F^x^i (left graphs) and axial- 
vector mesons Fa^ai (right graphs), for i = 1 (elastic case, solid line), i = 2 (dashed line), 
2 = 3 (dot-dashed line). 



Observe that, as — 0, the non-elastic form factors go to zero, while the elastic form factors 
approach one in this limit. As we increase g^, the form factors decrease at different rates. At 
small g^, F^i dominates, then F^i„2 takes over, then F„i„3, and so forth. The same situation 
occurs for the axial- vector case. Physically, this can be explained by the fact that, as the 
momentum transfer increases, the amplitude for producing heavier final states increases as 
well. 
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Regarding the large dependence of these non-elastic form factors, we can make an expan- 
sion similar to eq. (15.241) but for different initial and final states. If the following supercon- 
vergence relations hold, 



gy'^Qv^v^v^ — , gv^gy^a-^gi — 0. (5.36) 



n=l n=l 



we can again expect to find a g ^ behavior as in the elastic case. We calculated these sums 
from n = 1, 6, for the relevant states, finding: 

6 6 

^gvr^gv^v^v^ ~ -0.000105(GeV)^ ^^^n^^„„i„2 ^ -0.00316(GeV)^ 

n=l n=l 
6 6 

Y^gv^^gv-v^v^ - -0.00173(GeV)^ ^(7,n(7,„,i,3 ^ 0.0239(GeV)^ 

n=l n=l 
6 6 

'^g^ng^n^i^A 0.0155(GeV)^ y^^gyr^gv^a^a't ^ -0.185(GeV)^ . (5.37) 



n=l n=l 



These results indicate the validity of relations (I5.36p . The numerical errors in F^i^i and -Fai^i 
increase as we increase i. In order to check the large behavior, the transition form factors 
multiplied by are plotted in figure [9l along with the corresponding elastic form factors 
for comparison. We conclude that the transition form factors asymptotically approach the 
expected dependence, within the numerical errors. 



6. Conclusions 



In this paper we studied the spectrum of vector and axial-vector mesons, along with the 
pion and meson form factors, in the holographic model of chiral symmetry breaking of 
Kuperstein and Sonnenschein. Moreover, we found that vector meson dominance is a feature 
of electromagnetic interactions of mesons in the KS model. We compared our numerical 
results to the Sakai-Sugimoto model as well as expectations from QCD. Our results for 
the form factors are generally in good qualitative agreement with QCD results. The pion 
form factor obtained in our model agrees well even quantitatively with the data provided 
by the Jefferson Lab Fpi2 collaboration. This agreement is somewhat surprising, since our 
model is not dual to (large Nc) QCD, but rather a quiver gauge theory with gauge group 
SU{Nc) X SU{Nc). The meson physics is governed by open string excitations on the D7 
flavor branes and depends on the details of the D7-brane embedding. This embedding 
is in many ways similar to the embedding of flavor branes in the Sakai-Sugimoto model. 
However, the asymptotic behavior of the geometries is quite different. The asymptotic 
conformal invariance of the Kuperstein-Sonnenschein model may explain why our results at 
large are close to QCD. Another notable difference between the Sakai-Sugimoto model 
and the KS model, concerning numerics, is that, in many cases, convergence of certain 
integrals and summations is much slower and harder to achieve numerically, due to the slow 
fall-off of the KS wave functions for large z. Nevertheless, our results, especially for the 



23 



pion form factor, clearly demonstrate the viability and phenomenological potential of the 
Kuperstein-Sonnenschein model. It would be very interesting to improve our understanding 
of the differences between the D-brane embeddings with respect to QCD phenomenology. 
For example, one could try to classify possible D-brane geometries and the resulting five- 
dimensional effective actions and investigate their usefulness for phenomenological questions. 
We leave this question to future research. 

In a more recent paper, Dymarsky, Kuperstein and Sonnenschein (DKS) [30] study a similar, 
yet more complicated, model of chiral symmetry breaking. This model is based on the 
Klebanov-Strassler background [31] which can be considered a non-conformal generalization 
of the Klebanov-Witten background. It features confinement at low energies and flows (via 
a Seiberg duality cascade) to a theory close to M = 1 super Yang-Mills in the IR. With the 
addition of fundamental flavor degrees of freedom this becomes a QCD-like theory. We are 
presently investigating its phenomenology which we will report in a follow-up paper [32] . 
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